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STRUCTURE FUNCTION AND VARIABILITY MECHANISM OF QUASARS FROM SDSS STRIPE 82 
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Draft version 

ABSTRACT 

Theoretical predictions for the ensemble quasar structure function are tested using multi-epoch 
observations of Stripe 82 collected by the Sloan Digital Sky Survey. We reanalyze the entire available 
volume of the (7-band imaging data using difference image photometry and build high quality light 
curves for 7562 spectroscopically confirmed quasars. Our structure function includes ~ 4.8 x 10^ pairs 
of measurements and covers a wide range of time lags between 3 days and 6.9 years in the quasar 
rest frame. A broken power-law fit to this data shows the presence of two slopes ai = 0.33 and 
a2 = 0.79 with the break at ~ 42 days. The structure function compiled using only flux increases is 
slightly lower than that for variations of the opposite sign, revealing a slight asymmetry between the 
leading and trailing edge of a typical flare. The reality of these features is confirmed with monte-carlo 
simulations. We give simple interpretation of the results in the frames of existing theoretical models. 



Subject headings: quasars 



1. INTRODUCTION 



Several decades after the discovery of quasars, the 
physical origin of their variability is still not fully un- 
derstood. An unambiguous explanation of the intrinsic 
variations in quasar emission will likely provide the key 
to understanding the energy supply of quasars. The most 
frequently invoked model connects variability of quasars 
to instabilities arising in the a ccretion disk around a 
super-massive black hole (e.g. iKawaguchi et al.) 119981 : 
[Perevra ct al. 200(|). In the starburst model quasar 
light curves are a superposition of very frequent super- 
nova explosions in the host galaxv (e.g. iTerlevich et al.) 
Il992t ICid Fernandes et al.l Il99l lAretxaga et~all 119971 ). 
Gravitational microlensing by compact bodies within 
the host — such as normal stars or certain candidates 
for dark matter — also produces observable flux changes 
(e.g. Hawkins 2002 . 2007) . but their contribution to the 
observed level of variability remains uncertain. The mod- 
els can be distinguished based on their predictions for the 
dependence of the slope and the amplitude of the struc- 
ture function on the time-scale of variability in the quasar 
rest frame. Recently, individual quasar light curves have 
been successfull y modeled as a damped random walk 
(DRW) process (JKellv eraIll20Q9llKozlowski et aIll2Q10l 
iMcLeod et"aIll2010D . While the method allows one to 
study the structure function of individual objects as a 
function of black hol e mass, quasar lumi nosity, wave- 
length of observation (jMcLeod et al.l [20101 ) . it also fixes 
the slope of the power spectrum in the high frequency 
regime at —2 for every quasar. Moreover, the DRW 
model excludes the possibility that the structure func- 
tion of positive flux variations differs from the one for 
negative changes. In this paper we rely on more tradi- 
tional statistics using ensemble structure functions. 

A recent influx of survey data on variability of tens 
of thousands of QSOs is fueling numerous observational 
studies and calls for a new look at the proposed sce- 
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narios. The r nost recent paper s favor the disk in- 
stability model (iBauer et al.l l2009l Ide Vries et al.ll2005l 
Vand en Berk et aLll2004( l with the exception of lHawkinsI 
(,2002) who find evidence in support of the microlensing 
model. Here, we re-analyze the repeated ^-band imag- 
ing of Stripe 82 collected by the SDSS using difference 
image photometry and build an ensemble structure func- 
tion over a wide range of time-scales in order to test the 
theoretical predictions. Our results show that in their 
present form none of the considered models explains the 
behavior of the structure function on all time-scales. In- 
stead, a hybrid model including both a starburst and disk 
instabilities naturally follows from the data. 

The paper is organized as follows. In Section [2] we 
describe the data and a method of extracting light curves. 
Basic properties of the structure function are given in 
Section [31 In Section [H we build the structure function 
for the ensemble of quasars. We discuss the results in 
Section [5] and a summary is given in Section [S] 

2. DATA ANALYSIS 
2.1. Multi- epoch SDSS Imaging in Stripe 82 

For the purpose of this study we reanalyze all available 
g-band drift-scan images of the SDSS S tripe 82 in the fi- 
nal S DSS data release, hereafter DR7 (jAbazaiian et al.l 
I2009D . Quasar light curves were obtained us i ng difference 
image photometry ([Alard fc LuptonI Il998t lAlard[ l2000t) 
and a slightly modified version of the p hotometric soft - 
ware employed by the OGLE project ([Wozniakl |2000( ) . 
The original motivation for difference imaging approach 
was to f acilitate a sensitive search for strongly lensed 
quasars ([Kochanek et all 120061 : iLacki et al.ll2009t ). But 
the application of the method also returns high quality 
relative light curves of all sources in a given data set, 
which we use to study the variability of known quasars. 

The full database of Stripe 82 contains 303 runs cov- 
ering nearly 300 deg^ of sky along the celestial equator 
(-60° < RA < 60° and -1.25° < Decl < 1.25°). The 
area was scanned repeatedly during 10 years with some 
regio ns observed more than 80 times ([Abazajian et al.l 
l2C)09t) . The SDSS imaging system consists of 30 CCD 
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Fig. 1. — Spatial and temporal coverage of the SDSS Stripe 82 
imaging database. North runs are shown in blue and South runs 
in red. Note a much better time sampling in the later half of the 
survey. 



detectors organized in 6 raws (camcols), where a sin- 
gle CCD in each row is taking data in one of the 5 
filters ugriz. Since there are gaps between camcols, 
the width of the major stripe is split into North and 
South strips covered by independent scans. The SDSS 
database of Stripe 82 contains 162 North and 141 South 
runs. The spatial and temporal coverage of all runs is 
shown in Figure [1] where the blue and red lines corre- 
spond to the North and South runs. Each run is divided 
into 2048 x 1489 pix sections called fpC frames with a 
2049 X 128 pix margin on both ends ( Stoughton et al.l 
|2002[ ). The locations of the beginning and the end of 
each scan vary from one observation to another. This 
arrangement is not compatible with subtracting multiple 
images of the same area, so the runs must be "glued" 
back together and then subdivided using consistent field 
boundaries. The rebinning starts from RAmin, the first 
location covered by at least 20 runs (determined sepa- 
rately for North and South runs) . The rearranged frames 
of 2048 X 1489 pix are written as FITS files and the 
process continues untill we reach RAmax, beyond which 
fewer than 20 runs are available. 

Some of the data in Stripe 82 were obtained in non- 
photometric conditions and in relatively poor seeing 
(DR7). The performance of the difference image anal- 
ysis method is sensitive to the quality of the images, and 
therefore our next step is to reject images with poor see- 
ing, high background, and significant background gradi- 
ents. A small number of images taken toward the ends of 
the run have strong background gradients due to the Sun 
or Moon rising or setting. Occasionally, a problem with 
the bias level subtraction caused a jump in the count 
level between the left and right parts of the image corre- 
sponding to individual amplifiers on the CCD chip. Such 
images were removed from the analysis. A good quality 
reference image was prepared for each field using 10 im- 
ages with the best overall seeing, background, and trans- 
parency. This is done with the help of the SExtractor 
program (jBertin fc Arnoutal 19961 ) by requiring a median 



FWHM < 1.5" for point sources (CLASS_STAR > 0.7) 
and rejecting background maps with obvious problems. 
The fields with fewer than 20 remaining images were dis- 
carded. 

2.2. Difference Image Photometry of SDSS Quasars 

Interpolation of images to the same pixel grid was ac- 
complished using approximately 50 stars identified on all 
images of the field. The time base line of the Stripe 
82 data set approaches a decade. We found that a sig- 
nificant fraction of stars shows a detectable proper mo- 
tion producing time- variable residual s in the sha.pe of a 
dipole in the difference images (cf. iLacki et al.l 120091 ). 
In order to limit the impact of such problems on the 
quality of the final photometry, we selected the reference 
stars with the proper motion fi < 100 mas yr~^ from the 
iBramich et al.l (|2008( ) catalog in the overlap area (about 
50%). The spatial variability of the PSF matching ker- 
nel was computed using a third-order polynomial fit to 
model coefficients and a second-order polynomial for the 
differe ntial background Aga in, slow stars from the cat- 
alog of IBramich et al.l (|2008D were preferentially used to 
compute the solution. The DIA light curve is computed 
by adding difference fluxes measured with the DIA pack- 
age to the flux in the reference image. Object magni- 
tudes were converted to standard SDSS g magnitudes by 
adding a median offset between stars identified on the 
reference image and corresponding objects in the SDSS 
catalog. 

The photometric quality is demonstrated in Figure [2] 
showing the magnitude scatter versus baseline magnitude 
for a typical field. Here, the baseline is the median mag- 
nitude and the variability of the source is calculated as 
the 68-th percentile of absolute deviations from the me- 
dian. Black points correspond to the sources identified 
in that field and blue line shows median of the variability 
of the sources from all fields. 

In most fields the frame-to-frame accuracy reaches 7 
mmag for bright unsaturated sources. For a significant 
fraction of stars the photometric scatter is worse than 
that and limited by the residuals due to proper motions. 
The photometric uncertainty is 1% or better for sources 
with mg < 18 and around 10% for sources with rUg « 22. 
There are 52 quasars with nig > 22. Those were excluded 
from the structure function analysis Section |4|) due to 
large uncertainties in their light curves. In the next step, 
we match spectr oscopically confirmed q uasars from DR7 
Quasar Catalog (jSchneider et al.ll2010f ) with our photo- 
metric database. There are 9519 DR7 catalog quasars 
within the area covered by Stripe 82, and we can build 
DIA light curves for 7614 of them. There are several rea- 
sons why we miss a significant number of quasars. After 
exclusion of poor quality images many fields are left with 
fewer than 20 epochs and are removed from further anal- 
ysis. An average 10% of the image area is lost due to 
masking of very bright stars. Some fields are rejected for 
lack of stars that can be used to calibrate light curves. 
Finally, our pipeline occasionally fails on a small fraction 
of images with problems in fpC frame headers. 

An example quasar light curve resulting from this pro- 
cedure is shown in the Figure |3l There are 43 observa- 
tions in the light curve and time counts from the first 
epoch in Stripe 82 (run 94, fpC frame 12). The number 
of epochs per quasar in our sample varies between 20 and 
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Fig. 2. — Variability diagram in the g band for a typical field in 
SDSS Stripe 82. The magnitude scatter of all measurements for a 
given source is plotted against the baseline magnitude (median). 
The photometry is based on our reanalysis of the SDSS scans us- 
ing difference image photometry. Black points correspond to the 
sources identified by SExtaractor. Sources above the magenta line 
are considered to be variable. Blue line shows median variability 
of the sources from all fields. 
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Fig. 3. — Example difference imaging light curve for QSO 
SDSS 223154.15 + 004747.5 at z = 0.8409. Time counts from the 
first observation in Stripe 82. 



56 with the mean of 36. 

3. STRUCTURE FUNCTION 

The time samphng of our quasar hght curves is rela- 
tively sparse (see Figure [3]) and insufficient for the use of 
spectral methods such as the frequency-power spectrum 
analysis. Following the common practice in AGN vari- 
ability studies, we focus on the information contained in 
the structure function. A number of structure function 
definitions are in use differing in relatively minor details 
that are not important for this study. We adopt a struc- 
ture function (hereafter SF) in the form: 



S{t) 



Nir) 






(1) 



where Xi — x(ti) some random quantity at moment 
ti, the sum is taken over all measurements for which 
tj—ti = T and N{t) is the number of such data points. In 
short, the structure function is a time dependent stan- 
dard deviation of some random process. SF of white 
noise is flat and equal to v2ct^. Therefore, the SF of a 
random process is also flat on time-scales, where the am- 
plitude of the intrinsic variability is comparable with the 
measurement errors, as well as on the longest timescales 
where the r.m.s. stops changing as the length of the time 
interval keeps increasing. In the context of quasar vari- 
ability the random quantity is the magnitude (essentially 
the relative luminosity) rrii measured at time ti and r is 
taken in the quasar rest frame. For a given quasar we 
can write: 



S'ir) = ((m, - m,)2), ^ ^{{Sf, - 5J,f) 



1 



2ai 



{{5U-5L,Y)^^^, 



L2 



L2 



(2) 



where we converted to flux units. The flux consists of a 
stationary component / and a time dependent variable 
part 5fi, and Ur is the r.m.s. calculated for points sep- 
arated by time interval r. From here it is obvious that 
for two quasars at the same redshift and with equal a^ 
the brightest one will have smaller variability amplitude, 
S{t)^ at all time scales. 

The number of time intervals that can be constructed 
from an average of 36 epochs for a single quasar is 
36 X (36 — l)/2 = 630 and results in a very noisy struc- 
ture function"^. Therefore, we focus our analysis on the 
ensemble quasar SF, i.e. we assume that the ensemble 
variability reflects on average the variability of a single 
quasar. 

4. THE ENSEMBLE STRUCTURE FUNCTION OF 
QUASARS 

The final sample contains 7562 light curves. All pos- 
sible pairs of measurements within each light curve are 
collected for a total of 4,796,468 magnitude differences. 
After reducing the time lag to the quasar rest frame us- 
ing spectroscopic redshift s from the SDSS Quasar Cata- 
log ([Schneider et al.ll2010i' ). the data are binned in time. 
The bin size was adjusted to reach 20,000 measurements 
per bin. The resulting ensemble SF consists of 239 bins 
and is shown in Figure 2] (blue points). Error bars 
were estimated by dividing each timescale bin into 20 
smaller bins, evaluating the mean SF in each, and cal- 
cula ting the r.m.s. scatter of the resulting values (see 
also lde Vrieset al.l (|2005f ll. 

In order to verify that we detect quasar variability and 
guard against the possibility of contamination by unre- 
lated factors, we repeated this calculation for constant 
stars. For that purpose we select stars below the ma- 
genta line in Figure [5] with nig < 22. The result is shown 

^ Structure functions of individual sourc es can be effectively 
used for identification of quasar s (see e.g. ISchmidt et al.l 120101 : 
IPalangue-Delabrouille et al.ll20Til) . 
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Fig. 4. — Top panel: Raw ensemble structure function of 7562 
quasars from SDSS Stripe 82 (blue) compared to the structure 
function for ~ 1000 non- variable stars (green). The stars included 
in this comparison fall below the magenta line in Figure[2]| is shown 
by green points. Bottom panel: Noise-corrected structure function 
for an ensemble of quasars (blue) and the best fit broken power-law 
model (red). The oscillations for time-scales above ~ 100 days are 
an artifact of the time sampling. The model offers only an approx- 
imate description of the structure function over the full extent of 
the data. 



in the top panel of Figure |3] (green points). As expected, 
the SF of non-variable stars is flat over the entire range 
of time-scales. We conclude that the observed variability 
of quasars is intrinsic and not caused by unrecognized 
measurement errors. The stellar structure function ex- 
ceeds the noise level for the quasar sample because the 
stars selected for this comparison are fainter on average 
than the quasars. 

The quasar SF in the top panel of Figure |4] displays 
several slopes and breaks. It starts flat at r < 2 days, 
then at 2 < r < 10 days gradually transitions to an ap- 
proximate power-law increase over 10 < r < 40 days, 
and finally breaks to a shallower slope above r ~ 100. 
The flat portion of the SF is clearly due to noise and 
the transition to the flrst power-law section is due to 
SFs of single quasars with different mean brightness mg 
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Fig. 5. — Noise-corrected SF^ (blue) and SF^ ( green). Black 
oints show the g-band structure function from iWilhite et al.l 



II2008I) who used the SF^ estimator. Note a close agreement in 
both the shape and amplitude. 



reaching the level of noise at different r (see Equation[2]). 
The slopes seen at 10 < r < 40 days and t > 100 have a 
physical origin and are discussed further below. The time 
interval between 40 and 100 days is poorly sampled and 
difficult to interpret. This is a result of the SDSS scan- 
ning strategy for Stripe 82 that produces relatively few 
observations separated by 90-270 days in the observer 
frame (cf. Figure |3]) combined with the mean redshift of 
the sample ^1.5. 

The noise contributes to the SF equally at all 
timescales and must be subtracted in order to expose 
the intrinsic quasar variability. The noise-corrected SF 
is shown in the bottom panel of Figure |4] (blue points). 
The shape of the quasar SF is dominated by two slopes 
that we estimate by fitting a broken power law model to 
the data 



Sir) 



Pir/T^r 
/3(T/ro)"^ 



if T > To 
if T < To 



(3) 



We find ai = 0.33, a2 = 0.79, tq = 42.3 days, and 
P = 0.13 (solid red line). The reduced x^ of the fit 
is far from unity, i.e. Equation [3] does not provide a 
complete description of the data. However, the main 
goal here is to detect and estimate the two slopes in the 
overall shape of the SF for an ensemble of quasars, and 
compare those coefficients with theoretical predictions. 
We also estimate errors on the derived parameters using 
'jackknife' resampling. Namely, we randomly select 1/3 
part of the sample, drop it, and fit ensemble SF built 
from the rest. We repeat the procedure 100 times and 
estimate the error on the given parameter as r.m.s. of 
proper distribution. Such estimation gives very small 
errors (0.01, 0.02, 3.9 days, and 0.01 for ai, a2, tq, and 
/3 correspondingly) saying that the shape of ensemble SF 
is quite stable to small variations of the sample. 

For compa rison we plot the (?-band SF from 
IWilhite et al.l ([2008) (black points in the bottom panel of 
Figure |4|) . The overall agreement in shape between the 
SFs is very good and the discrepancy in no rmalization is 
explained by the different SF type used in IWilhite et al.l 
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Fig. 6. — Light curve simulations. A random simulated light 
curve (red) is compared with a real light curve of one of the quasars 
in the sample (blue). The simulated light curves were generated 
from a power spectrum in the form of a single power-law with the 
slope —1.66 (top) and —2.58 (bottom) corresponding to structure 
functions in Figure [7] 



(|2008f ). There are two types of SF commonly used in lit- 
erature. Following Bauer et al. (2009,) we refer to them 
as SF-^ and SF^: 



S^{T)^J^{\Am\y 



(4) 
(5) 



T he structure function in the form SF^ was introduced 
bv ldi Clemente et al.lll996l Our noise-corr ected SFs and 
the g-band SF^ from lWilhite et alj ()2008[) are shown in 
Figure [51 There is a nearly perfect agreement in both 
shape and normalization between the B-type SFs derived 
from the two data sets. 
The interpretation of results obtained using 

SF as a statistical tool is not straightforward. 

lEmmanoulopoulos et alJ ()2010[ ) show that spurious 
breaks can appear in the SF due to sparse time sampling 
and windowing effects. Moreover, the slopes derived 
from SF modeling are typically much more uncertain 
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Fig. 7. — Confirming the presence of two slopes in the quasar 
structure function. The simulated light curves were generated from 
a power spectrum in the form of a single power-law with the slope 
— 1.66 (blue) and —2.58 (green) corresponding to the upper and 
lower portions of the best fit model (red line) from Figure |4] 



than suggested by naive error estimates. Does it mean 
that the presence of two slopes in Figure |4] is an artifact? 
To answer this question and validate our results we 
perform a set of monte-carlo simulations. The goal of 
this calculation is to verify that for a set of light curves 
generated from a single power-law spectrum with a given 
value of a and having the time sampling of real quasars, 
the algorithm of computing SF does not introduce any 
significant breaks or biases in the estimated slope. 

The simulation is b ased on the algorithm of 
iTimmer fc Koenigl (|1995l ) and reproduces both the cor- 
rect power spectrum as well as the phase mixture. The 
power spectrum in the form P{f) oc j-(2a+i) corre- 
sponds to SF S{t) cx t". We simula te light curves 
accord ing to a prescription in lEmmano ulopoulos et all 
(|2010D . The simulated time series are then normalized 
to have the same r.m.s. and mean as the real light curves. 
Examples are shown in the Figure |6l 

We simulate two sets of light curves with the power-law 
spectral density (—1.66 and —2.58), select points on the 
simulated light curves corresponding to the time sam- 
pling of a real light curve, and apply proper rest frame 
correction. In Figure [7] the simulated SFs are compared 
to the best fit model for real quasars. The slopes of 
the simulated SFs do not exactly match the slope of the 
model from which they were generated, i.e. the time 
sampling slightly biases the measured shape of the SF. 
However, none of the simulated SFs is even close to a bro- 
ken power law with ai — 0.33 and a2 = 0.79. Therefore, 
the large observed difference in the SF slope between the 
short and long timescales is not due to the time sampling, 
but rather of a physical origin. Possible interpretations 
are discussed in the next section. 

The ensemble SF considered here is built from quasars 
with different intrinsic characteristics such as black hole 
mass or luminosity. Therefore, it is interesting to check 
whether two slopes in SF are characteristic for the hole 
range of, say, black hole masses or have dependence on 
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TABLE 1 
Best fit parameters 



Subsample 



Ql 



«2 



TO, days 



/3 



low mass 0.27 ± 0.02 0.72 ± 0.04 55.0 ± 16.7 0.15 ± 0.02 
mean mass 0.30 ± 0.01 0.74 ± 0.03 57.3 ± 6.6 0.15 ±0.02 
high mass 0.37 ± 0.02 0.94 ± 0.09 36.6 ± 15.6 0.10 ± 0.02 



mass. Using black hole mass estimations from lShen et al.l 
(|2008f ) , we select three subsamples of quasars with black 
hole masses in ranges < 5 x 10^ Mq (1953 sources, low 
mass range), 5 x 10* -^ 1.5 x 10^ Mq (2108 sources, mean 
mass range), > 1.5 x IO^Mq (1276 sources, high mass 
range) and build SF for each range of masses. All three 
SFs show existence of two slopes (see Figure 13 top 
panel). We fit SFs by broken power low model (Equa- 
tion [5]). The best fit parameters and formal estimation 
of errors by 'jackknife' resampling are given in the Ta- 
ble[T] In order to have better understanding of the uncer- 
tainties in the parameters we also plot regions showing 
minimal and maximal values of the fits for all three SFs 
(see Figure HI bottom panel) . As it is seen from the Fig- 
ure [SJ the SFs slopes tend to increase with the growth of 
the mean subsample mass but still have common values. 
More subtle study is needed in order to show that the 
trend in slopes as a function of black hole mass is real. 
Here we only conclude that the presence of two slopes in 
the ensemble SF is independent on the considered black 
hole mass range. 

5. DISCUSSION 

We consider three processes frequently invoked to ex- 
plain the variability of quasars: instabilities in the ac- 
cretion flow around a supermassive black hole, super- 
nova explosions related to the starburst phenomenon, 
and gravitat ional microlensing by compact bodies in the 
host galaxy (jKawaguchi et alj|l998HHawkin!Jl2002f) . SFs 
of light curves predicted by each model are character- 
ized by different slopes. The expected SF slopes are 
0.25 ± 0.03, 0.44 ± 0.03, 0.83 ± 0.08 respectively for vari- 
ability generated by m iscrolensing, dis k instability, and 
starburst mechanisms (jHawkinj l2002t ) . The SF slopes 
found in Section |4] are 0.79 and 0.33 correspondingly for 
time lags below and above 42 days. Taken at face value, 
the first slope is consistent with the starburst model, and 
the second slope falls half way between the predictions 
of the disk instability model and the microlensing model. 
Note that SF slopes in the range 0.30-0.35 have been ob- 
tained in other works, but the autho rs still attribute d the 
variability to disk inst abilities. iBauer et al.l (|2009l ) and 
iMeusinger et al.l (|2011[ ) present an overview of recent re- 
sults. 

Asymmetries in the SF provide an additional test of 
model predictions. The SF that only includes those 
pairs of measurements for which the fiux increases with 
time, 5*+, may be differe nt from S-, the SF cha racter- 
izing a decrease in flux. iKawaguchi et al.l ()1998l ) show 
that in starburst models S+ > S-. Disk instabilities 
produce 5_ > 5+ and microlensing is symmetric, i.e. 
iS*- = S^ when a veraged over sufficiently long time in- 
tervals (jHawkinsI I2OO2.) . The difference between S- and 
5*+ depends on the parameters of the model and there- 
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Fig. 8. — Top panel: Ensemble SF for quasars with black hole 
masses lying in low mass range (red points), mean mass range 
(green points ), and high mass range (blue points) (see text for 
definition of ranges). Thin black lines show best fits for each SF. 
Bottom panel: Regions showing the scatter of the best fits obtained 
during 'jackknife' resampling. Red lines show the scatter of the low 
mass range SFs best fits, green region and blue lines show the same 
for mean mass and high mass ranges correspondingly. 



fore S- ~ 5*+ is still possi ble in both the starbur st and 
the disk instability model (jKawaguchi et al.lll998D . 

Noise-corrected S+ and S- functions are shown in Fig- 
ure [9] with blue and green lines respectively. For long 
time lags in the range 300-1600 days, where the SF slope 
is ~ 0.33, we have S- > S+, as predicted by the disk 
instability model. The significance of the difference is 
not very high, but allows us to exclude the microlensing 
model. The difference disappears for time-scales below 
100 days, where the slope is consistent with the starburst 
model. This may suggest that a typical supernova rate 
in a starbursts in the host galaxies of quasars is as high 
as ~ 100 yr~^ and effect ively washes out any de tectable 
asymmetry (see Fig. 6 in lKawaguchi et al.l (|1998[ )V How- 
ever, the rate should be treated with caution since it is 
not derived from direct measurements. Moreover, the 
predictions for 5+ and S- were derived from simulations 
of a single variability mechanism. The addition of an- 
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Fig. 9. — Testing the asymmetry in the quasar structure function 
between the leading and traihng branch of a typical flare. Flux 
increases contribute to S-^. (blue) and flux drops contribute to S— 
(green). 



other variability process can wash out the difference be- 
tween 5*+ and S-. This is an interesting topic for future 
research. The difference in S+ and S- was also investi- 
gated by Ide Vries et all (|2005l) and iBauer et all (|2009t ). 
The former study found tentative evidence that 5*+ > S- 
(light curves have fast rise and slow decay as in starburst 
model) for 400 days < r < 1500 days, while in the lat- 
ter work S+ ~ 5- over the full range of measurements 
(r < 1000 days). The explanation of the difference in 
results is difficult because each paper is based on en- 
tirely different data and requires further research. The 
Ide Vries et al.l ()2005D analysis relies on cross survey flux 
comparisons. 

The introduction of the supernova variability mech- 
anism may be redundant as the shape of the SF can 
be explained by detailed variability mechanisms act- 
ing in accretion disks. For instance, the PSD of X- 
ray light cu rves of some a. c tive g alactic nuclei have 
two (e. s. [Summons et alj (120071 )) or more distinct 
slopes (jMcHardv et al.ll2007[ ). Considering that the op- 
tical emission tends to follow the X-ray flux with some 
time delay one may expect that PSDs of optical light 
curves will display two slopes. Moreover, the SF shape 
found here is in qualitative agr eement with the model 
proposed in lArevalo et al.l (|2009[ ). where large amplitude 
optical variations on the time scales of hundreds of days 
are attributed to the fluctuations in accretion rate and 
short-term (order of days), small amplitude variations 
arise due to reprocessing of X-rays. 

6. SUMMARY 

We applied the image differencing technique to the en- 
tire g-band imaging data set of SDSS Stripe 82 and con- 
structed high quality light curves for 7562 spectroscopi- 
cally identified quasars with no less then 20 epochs. The 
variability analysis was performed using the ensemble 
structure function. While the shape and normalization of 
the quasar structure function derived here are in perfect 
agreement with the results of lWilhite et al.l (|2008D based 
on an earlier SDSS data release, our structure function 



is estimated from a much larger data set and covers a 
substantially wider range of time lags from 8 hours to 
6.9 years in the quasar rest frame. 

We find that the ensemble SF reveals two distinct 
power-law slopes with the break around 42 days and 
confirm the presence of these features with monte-carlo 
simulations. The presence of two slopes in the ensem- 
ble SF is independent on the black hole mass range of 
quasars used to build it. The slopes estimated from 
a broken pawer-law fit to the data are ai ~ 0.79 for 
r < 100 days and a2 — 0.33 for r > 300 days. Using 
predictions from the theoretical models, the variability 
of quasars on times scales r < 100 days can be explained 
by a starburst model. The slope of the structure func- 
tion on time-scales longer than about 3 months is rather 
close to models where variability is explained by gravi- 
tational microlensing. However, on time-scales r > 300 
days we detect a significant asymmetry in the SF charac- 
teristic for disk instabilities that rules out a substantial 
contribution of microlensing. It is also possible to explain 
the shape of SF without invoking the starburst model, 
but using the model, where short term optical variations 
are caused by X-rays reprocessing and long therm varia- 
tions originate from optic al emitting regions of accretion 
disk (jArevaloet al.ll2009t) . 

The existing SDSS Stripe 82 data offer a great oppor- 
tunity to better understand the origin of variability in 
quasars. The main limiting factor at present is the lack 
of detailed theoretical predictions that go beyond the ba- 
sic characteristics such as the SF slope and include the 
detailed shape and normalization of the SF, as well as 
asymmetries and color information. 
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